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Abstract 

We discuss "mode-entangled states" based on the optical transverse modes of the optical field 
propagating in multi-mode waveguides, which are classical simualtion of the quantum entangled 
states. The simualtion is discussed in detail, including the violation of the Bell inequality and the 
correlation properties of optical pulses' group delays. The research on this simulation may be im- 
portant, for it not only provides useful insights into fundamental features of quantum entanglement, 
but also yields new insights into quantum computation and quantum communication. 
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Introduction 



So fai\ there is interest in research on classical wave analogs of the Schrodinger wave 
function QDD . It is well known that, in the paraxial approximation, the transverse modes 
of an optical field obey a propagation equation which is formally identical to the Schrodinger 

n 

equation with the time replaced by the axial coordinate The transverse modes of the 
optical field propagating in a waveguide with a parabolic refractive index profile are formally 
identical to quantum harmonic oscillator wave functions. Some efforts have gone into re- 
searching on classical wave analogs of quantum mechanics, including analogs of Fock states 
and measurement of Wigner phase-space distributions for classical optical fields which can 
exhib.t negative regions^. However, research on epical analogs has been lifted 
principally to measurement of first order coherence, i.e., single-particle states. Classical- wave 
ana>ogs of high order coherence (quantum entanglement), i.e., mu.tipart.de states, have been 
seldom studied |7j. The quantum entanglement, which describes nonlocal quantum correla- 
tion between different degrees of freedom especially separated particles, is regarded as the 
inherent feature of quantum theory ^]. The quantum correlation has been shown in the 
correlation measurement of the entangled state, and a criterion has been given by the viola- 
tion of the Bell inequality In recent research, the quantum entanglement is considered 
as a key property to realize the quantum computation 1(| and quantum teleportation 
which makes the quantum entanglement strongly attracted to researchers. 

In this paper, we will propose "mode-entangled states" based on the optical transverse 
modes of the optical field propagating in multimode waveguides. It is well-known that the 
quantum entanglement is the characteristic of the quantum theory with no classical analog. 
Therefore, the "mode-entangled states" should be interpreted as the classical simulation 
of quantum entanglement using the optical transverse modes of the optical fields. The 
classical simulation will be discussed in detail, including the violation of the Bell inequality 
and the correlation properties of optical pulses' group delays. In [3], a full optical scheme 
to perform quantum computation is proposed, based on the optical transverse modes in 
multimode waveguides. The proposed C-NOT gate has the potential of being easily realized 
since it is based on optical waveguide technology and can be constructed by using Mach- 
Zehnder interferometer having semiconductor optical amplifiers (SOAs) in its arms. The 
SOA can provide a very large Kerr-like nonlinearity even at relatively low light intensities 
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and can avoid the intensity attenuation excited by two-photon absorption 15]. Therefore 
SOAs have been used extensively as nonlinear elements in optical switching and wavelength- 
conversion devices [l^. But whether the scheme would be capable of implementing the 
quantum computation depends on whether the C-NOT gate proposed in |l2j can generate 
the "mode-entangled states". Given all these, the research on the similarity between the 
"mode-entangled states" and the quantum entangled states may be important, for it not 
only provides a useful insight into fundamental features of quantum entanglement, but also 
yields new insights into quantum computation and quantum communication. 

The paper is organized as follows: In Section HI we will discuss the analogies between 
optical transverse modes in a multimode waveguide and quantum Fock states. In Section 
HTl the superposition of transverse modes in a random waveguide is analyzed. In Section 
IIII1 the Bell inequality as a criterion of the existence of mode-entangled states is deduced. 
In Section IIV1 an analysis of the mode-entangled states in random waveguides and the 
correlation properties of group delays is discussed. Finally, we summarize our conclusions 
in Section El 

I. ANALOGS OF QUANTUM FOCK STATES USING OPTICAL TRANSVERSE 
MODES 

Considering a weakly guiding, symmetric slab waveguide, an optical field in the prop- 
agation direction, longitudinal z direction, is restricted within the core region, which has 
the higher refractive index (RI) compared with that of the cladding. By using the Fock- 
Leontovich paraxial approximation, the Maxwell equations for the monochromatic electric 
field component can be reduced to the equivalent Schrodinger equation. The reduced field 

^ (x, y, z) = y/n^E (x, y, z) exp (-ik n dz^j (1) 
satisfies the following equation 

1% = ( V ** + + \ K - n 2 (*, y, z)] * = m (2) 

where E (x, y, z) is the monochromatic electric field component, the geometry of the waveg- 
uide is defined by the RI profile n(x,y,z), n = n(0,0,z), k = 2n/\ (with A being the 
wavelength in free space) and £ = j^dz/riQ. Therefore, when the RI profile n(x,y,z) is 



parabolic, the equation (J2J) is similar to the Schrodinger equation for a quantum harmonic 
oscillator. Here, instead of Plank's constant h, we have the vacuum light wavelength A. And 
the variable £ plays the role of time. Following [l^j], coordinate and momentum operators X 



and P 



(d/dxi 



1, 2 denote x and y respectively) can be introduced. These oper- 



ators obey the standard commutation relations 



X,, Pj 



(i/k)5ij and uncertain relations 



(AX,) 2 (AP,) 2 > (1/4P), where (AX) 2 = (xf) - (x) and (AP) 2 = (P l 



Pi 

When the optical field propagates in the waveguide that is ^-invariant, in other words, 
the RI profile is uniform along z, the propagation can be equivalently described in the 
time-independent Schrodinger equation 



H^ n (x,y) = uj n \V n (x,y) 



(3) 



where ty n (x, y) are a set of eigenmodes corresponding to a set of discrete eigenvalues uj n . It 
leads to the expression of the monochromatic electric field component E (x, y, z) as follows 

E(x,y,z) = J2C n e-^ n (x,y) (4) 

n 

where the propagation constants f3 n are given as 

f3 n = k(nl-2u n ) 1/2 (5) 

As described in [3], such eigenmodes \V n (x, y) are similar to the quantum Fock states. Here 
we introduce the annihilation operators fij 
at = 



k/2 (X + iPi) and the creation operators 



k/2 (Jti — iP^j that obey the boson commutation relations 



Oj , (Ij 



fail Oj] 



0. Application of the creation and annihilation operators to the Fock states yield 

a + \n) = y/n + 1 \n + 1) 
a \n) = y/n \n — 1) (6) 
a + a \n) = n \n) 

where |n) are the eigenmodes \V n (x), n = 0,1,2,..., (for simplicity only x direction is 
considered) . 

It is well known that random perturbations of the geometry of multimode optical waveg- 



uides cause fluctuations of the average arriva 
optical pulse propagating in the waveguides 



time (group 



18] 19] 



20] 21] 



delay) and spread (dispersion) of 



221 ] . In general, the information 



given for the description of optical field propagation in a random waveguide by means of the 



field ^ (x, y, z) in Eq. ((TJ) is not complete. The optical field propagation may be generally 
described by means of the density matrix formalism Q 



P = Y,Pmn\™) H (7) 

which satisfies the Liouville equation 



nm 



The density matrix possesses the usual properties: Trp = l,Trp 2 < 1 (the equality is true 
for pure states). The expectation value of any operator Q is given by the trace of the product 
of p and Q: (Q) = Tr (pQj . Therefore, the utilization of the density matrix formalism seems 
to be useful for describing a superposition of modes in the random waveguide. 



II. SUPERPOSITION OF TRANSVERSE MODES IN A RANDOM WAVEGUIDE 

In [3], a full optical method based on the transverse eigenmodes is proposed to perform 
the quantum computation, in which TE mode and TEi mode in dual-mode waveguide are 
used as qubits to represent logical and 1. In this section, we will use the density matrix 
formalism in the analysis of superposition of these two modes (TEo mode and TEi mode) 
in a random waveguide. 

The superposition of the modes can be described as 

* (x, y, z) = Coe^* |TE > + C x e~^ z |TEx) (9) 

where j3 Q and fi x are the propagation constants of the modes |TE ) and |TEj), respectively. 
In the dual-mode waveguide, the coupling of TEO mode and TEI mode is similar to a 
two-level system. Thus to describe this kind of coupling, we introduce the Hamiltonian 

H = O (a+a + \)+Pi (b + b C ab a + b + C* ab b + a (10) 

where a + and a are the creation and the annihilation operators of the mode |TE ), and b + 
and b are the creation and the annihilation operators of the mode |TEi). In the random 
waveguide, the random coupling among the guided modes will be caused by the perturbations 
in the waveguide geometry. Here we introduce a coupling coefficient to describe the random 
coupling. The coupling coefficients are functions of z coordinate that measures distance 
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along the waveguide axis. In random waveguides, the coupling coefficient assume the form 

C ab = K ab f(z) (11) 

where K ab is independent of z. The function f(z) often describes the actual shape of the 
deformed waveguide boundary or the bent waveguide axis. And it is supposed to be a 
stationary random variable whose correlation function is assumed to be Gaussian 

(/ (z) f (z - «)) = <?e-W (12) 

where (. . .) denotes the average over an ensemble of random realizations, a is the variance 
and D is the correlation length of f(z). 

To study effects on the superposition of transverse eigenmodes caused by this kind of 
randomicity, we introduce the density matrix p. If we rewrite the modes |TEo) and |TEi) 

as 





|TE ) =1*1* l TE i) 



V 



1 

the density matrix p of the mode superposition (jHJ) can also be rewritten 



(13) 



\C*' C C^ 
dC *e- iA ^ |Ci| 2 

where A/3 = f3 1 — (3 . The Liouville equation (|SJ) can be written in the equivalent form 

= {(3m _ Pmn {z) + {[V {z) , p {z)])mn (is) 

where m,n G {0, 1}, V (z) = C a ba + b + C* b b + a. By using the method mentioned in j^, after 
the optical field propagates a distance of L in the random waveguide, the density matrix p 
can be described as 



P 



where 



(l + e - 2 ^) &f + (l - e" 2 ^) C C* ie ^ + ^ L \ 



d^el- 4 ^-^ (l + e~ 2 ^ L ) J^L + (l - e 



(16) 



7 = ^De-^M 2 \K ab \ 2 (17) 



k = Im 



v^a 2 De-( DA W 2 erf (iDA/3/2) \K, 



ab | 
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We have assumed so far that A/3 ^> 7, Af3 ^> k and the waveguide is lossless. From Eq. 
f!16|) . we can anticipate if L — > 00, Tr p 2 — > |, which shows the evolution from the coherence 
(pure state) superposition to the incoherence (mixed state) superposition caused by the 
perturbation in the random waveguide. 

In order to distinguish between the coherence superposition and the incoherence super- 
position, we propose a scheme by analyzing the symmetries of the eigenmodes. In Fig. 1, 
the profiles of the modes |TE ) and |TEi) are shown, where |TE ) is symmetric and |TEi) is 
antisymmetric. Y-splitter is a device to split one light beam into two beams. If the pertur- 
bation of Y-splitter's transition is slight (i.e. the change of propagation constant A/3 ~ 0), 
the splitter can split the beam with extremely low power loss. When |TEo) and |TEi) are 
launched in the Y-splitter, |TEn) is split into two symmetric parts, while |TEi) is split into 



24j . Therefore, the output states of Y-splitter's two branches are 



two antisymmetric parts 
given explicitly by 

|+) = -L(|TE ) + |TE 1 )) (18) 

I") = -^(ITEoMTEx)) 

When the input of Y-splitter is given at a coherence superposition, \^fm) — 
{e~ ie |TE ) + e %e |TEi)V via expanding the input states in terms of the output states 
of the two branches of the splitter, we got \^i n ) = cos# |+) — isinO \ — ), from which the 
intensities of two branches can be obtained 

\(+\V m )\ 2 = cos 2 6 (19) 
|(-|M/ m )| 2 = sin 2 fl 

Note that the phase 9 of the input field \^i n ) will cause the variation of intensity in the 
Y-splitter's two branches. When the input of Y-splitter is given at an incoherence superpo- 
sition, the relationship of the phase between the modes |TEo) and |TEi) is uncertain. The 
density matrix of the incoherence superposition is given as 

p= ^ n |TE n )(TE n | (20) 

71=0,1 

where W n = \ are the probabilities for the two modes |TE ) and |TEi). Then the output 
average intensities in two branches are (+|p|+) = (— \p\— ) = \- It shows that when the 







1 e 2ie \ 
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1 -e 2ie 


2 




e~ 2i9 1 



input is at the incoherence superposition, no matter how the phase of the input changes, the 
output intensities will stay invariable. That is the essence of our scheme of measuring the 
intensity difference between Y-splitter's two branches to distinguish between the coherence 
superposition and the incoherence superposition. 

Now, we apply the Y-splitter to the analysis of the superposition state that is the evolution 
of a coherence superposition state after propagating in a random waveguide with distance 
of L. We define operators to represent the operations of a phase controller and Y-splitter 

I I 1 e 2W \ 

I + {6) = P + (0)|+)<+|P(0) = -[ _ j (21) 

t(9) = P+ (6) |-) (-| P (6) 

where P (6) denotes the phase control of the input state 

P(0)|tf) = P (9) (Co |TE ) + Ci |TEi)) (22) 
= C e- ld |TE ) + C x e ie |TEi) 

The control of the phase difference between the two modes |TE ) and |TEi) can be achieved 
by properly changing the RI of the core layer j^J. 

After the input at a mode superposition state propagates a distance of L in the random 
waveguide and passes through a phase controller and a Y-splitter, the intensity difference 
between the output waveguides can be obtained by using the density matrix p in Eq. (fTEJ) 

(/+ (6) - r (6)) = Tr {p [/+ (6) - r (6)} } (23) 

When C = d = ±, (/+ (6) - t (6)) = e~^ L cos [29 + (A/3 + k) L\. It is already obvious 
that the perturbation in the waveguide geometry may cause the evolution from the coher- 
ence superposition to the incoherence superposition, namely decoherence, which leads to 
disappearance of the intensity difference between two outputs of the Y-splitter. 

By using beam propagation method (BPM) we calculate numerically the behavior of 
the mode superposition in the dual- waveguide when it is propagating in the phase controller 
and the Y-splitter discussed above. The results are shown in Fig. 2, which show that the 
intensities of Y-splitter's two branches vary by changing the RI of the core layer. 
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III. BELL INEQUALITY OF OPTICAL TRANSVERSE MODE ENTANGLE- 
MENT 



As shown in |3[, in the Wigner distribution, the optical transverse modes are similar to 
quantum Fock states. However, such similarities are confined to first order coherence (such 
as single-particle). The higher order coherence (such as multiparticle) is regarded as the 
inherent feature of quantum phenomena, which is nonlocal quantum correlation shown in 
the correlation measurement of the quantum entangled state. The criterion of the existence 
of the quantum entanglement is given by the violation of the Bell inequality. In this section, 
we will discuss the correlation properties of the optical mode entanglement that is classical 
simulation of the quantum entanglement. 

We assume that a kind of "mode-entangled states" can be generated by means of some 
kind of interaction between the optical fields propagating in multimode waveguides (e.g. the 
C-NOT gate proposed in |3]). The mode-entangled states are given as 

$±) = (|TE ) C |TE ) t ± |TE 1 ) c |TE!> t ) (24) 

Vf) = ^ (|TE > jTE!> t ± |TE!> jTE ) t ) 

where c and t represent the control and the target fields, respectively. The states in each 
waveguide are a mode superposition, but they are different from a product state 

|* a > = \ (|TEo> c + ITE^J (|TE ), + |TE x ) t ) (25) 



The difference of 



$f) (ttf)) and |* 2 ) can be obtained not by measuring a single field, 



but by the correlation measurement of the control and the target fields. This correlation 
measurement can show that the two entangled fields are impartible to some extent. Similar 
to quantum entanglement, the violation of the Bell inequality is also used as the criterion 
of this impartibility. 

To perform the correlation measurement of the Bell inequality, a mode analyzer with a 
phase controller is required. The construction mentioned in the last section consisting of 
a phase controller and a Y-splitter meets the requirements. Following Eq. (|2ip. we define 
operators if and if to represent the mode analyzer's operations on the control and the 
target fields, respectively 

it - k = P + {Ox) (l+) c (+\ c - He HJ p (Oi) (26) 



2i0i 



p + (02)(\+) t (+\ t -\-) t n t )p(02) 



= e^ 2 |TE ) t (TEi| t + e~ Zl ° 2 |TE 1 ) i (TE | t 

where P (#1) and P {9 2 ) represent the phase controllers on the control and the target fields, 
respectively. As discussed in section |TTJ when 9\ and 9 2 change, the output intensities of the 
Y-splitters will vary correspondingly. 

Based on the correlation analysis, we propose an experimental scheme, shown in Fig. 3, 
in which the mode-entangled state is generated via the C-NOT gate proposed in The 
input of control field is given at the mode superposition 4^ (|TE ) + |TEi)) and the input 
of target field is given at the mode |TEo) or |TEi). Then the output fields of the C-NOT 
gate are sent to spatially separated mode analyzers represented by if and I 2 ■ The detected 
photocurrents of the mode analyzers' outputs are passively subtracted and monitored on 
a spectrum analyzer (SA) to check for correlations. Therefore, the correlation function is 
given by 

"it-iiMti-k 



E\ 



,02) 



(27) 



Substituting and |^ 2 ) m to Eq. (|2*7|). respectively, we obtain the correlation functions 
of the two states 



E. 



9 2 ) 



If ~ II 



n + n 



(28) 









(if - 12) 


$+} 








(if + k) 





Eq, 2 (9i, 62) 



cos (2#i + 20 2 



(if - II 



if + in [I2+1 



(29) 



(*2l 




(if - k) 


l*2> 


(*2l 


(If+Il) 


(if + 12) 


1*2) 



= cos (20i) cos {29 2 ) 
Then we substitute the correlation functions above into the Bell inequality 
\B\ = \E 9 2 ) -E(9 1 , 9' 2 ) + E (9[, 9' 2 ) + E (9[, 9 2 )\<2 
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(30) 



This particular Bell inequality is known as Clause-Horne-Shimony-Holt (CHSH) inequality 
0| . For the entangled state when we choose 6\ = n/8, Q\ = — 7r/8, 9 2 = 0, 2 = 



we get 1 5 1 = 2^2. Obviously, by proper choice of the phases 0\ and 62 in Eq. (|2l)|) . the 
correlation of the analyzers can exhibit a maximum violation of the Bell inequality \B\ > 2. 
However the violation never occurs for the product state l^)- 

We have simulated numerically the scheme shown in Fig. 3 by using BPM. The result is 
illustrated in Fig. 4. Due to the limitation of the simulation, we can't get the correlation of 
the control and the target fields. Therefore, the experiment is necessary to validate whether 
there are mode-entangled states similar to quantum entangled states. 



IV. OPTICAL MODE-ENTANGLED STATES IN RANDOM WAVEGUIDES 



As shown in section ing the Bell inequality will be violated in the correlation measurement 
of a mode-entangled state. However, we can see from Eq. (|'23j) that the violation will vanish 
due to perturbations of the random waveguides. The perturbations cause fluctuations of the 
average arrival time (group delay) and spread (dispersion) of optical pulse propagating in the 
random waveguides. In this section, we will further discuss the difference of the entangled 
states 



$1/(^1 /) and the product state ^2) by analyzing the correlation properties of the 
group delays. And this difference is another proof of the existence of the mode-entangled 
states mentioned in the last section. 

After the control and the target fields of the mode-entangled state $^ propagate respec- 
tively in two random waveguides with the same random characteristic (the variance a and 
the correlation length D) and the same distance of L, the density matrix p can be described 
as 



1 





Q Q e 2[i(A/3+ K )- 7 ]L \ 






V 



(31) 



/ 



e 2[-i(A/3+ K )- 7 ]L q q 1 

Due to the perturbations of the random waveguides, the coherent properties of mode su- 
perposition of the control and the target fields will decay exponentially with increasing the 
distance of L, until the whole state evolves to an incoherent superposition of |TE ) C |TE ) t 
and |TEi) c |TEi) r Similarly, the density matrix p of the product state ^2) propagating in 
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the random waveguides can be described as 
/ 



P* 2 



I e [i(A/3+K)- 7 ]i e [i(Ap+K)-i\L e 2[i(A/3+ K )- 7 ]L \ 

e [-i(A/3+K)-j}L i e -2yL e [i(Ap+ K )-y]L 

e [-i(A/3+«)-7]i e -2 7 L ! e [i(A/3+«)- 7 ]i 

g 2[-i(A/3+K)- 7 ]L e [-i(A/3+ K )- 7 ]L e [-i(A/3+«)- 7 ]£ ^ 



(32) 



When L — > oo, the state ^2) evolves to an incoherent superposition of |TE ) C |TE ) t , 
|TEi) c |TEi) t , |TEi) c |TE ) t and |TE ) C |TEi) t , which is obviously different from the evolu- 
tion of the mode-entangled state $1"). And the difference can be shown by a correlation 
measurement of group delays. 

The group delay of an optical pulse in a waveguide can be expressed as 



Ld{3 
c dk 



(33) 



where j3 is the propagation constant, c is the light velocity and L is the length of waveguide. 
If we introduce the group delay operator f whose eigenvalues are the group delay r, the 
average arrival time of the optical pulse can be obtained as follows 



(r(L)) = Tr(pf) 



(34) 



To study the correlation measurement of the group delays, we define the correlation function 
between the group delays of the control and the target fields as 



(Ten) = ((r c -(r c ))(r t -(r t ))) 
= (r c T t ) ~ (t c ) (f t ) 



(35) 



Substituting Eqs. (jSH) and (|3^|) into Eq. (|33|) . the correlation functions of the entangled state 
and the product state |^ 2 ) are obtained 



(Tc,T t )*+ = Tr (p $ +f c f t ) - Tr (p cl f c ) Tr {p tl r t ) 

fd^ _ dpA 2 
= 4? [~dk ~dk) 

= 7 [*i ( L ) ~ ^ (L)} 2 



(36) 



(Tc,T t )u a = Tr [p^r c f t ) - Tr (p c2 f c ) Tr (p t2 f t ) 
= 



(37) 
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where the reduced density matrices p cl , p tl , p c2 and p t2 are the partial traces Tr t (p$+J) 
Tr c Tr 4 (p# 2 J and Tr c (p^, 2 ), respectively. Here t (L) and ti (L) are the propagation 

time in the waveguide with distance of L for the modes |TEo) and |TEi) respectively. From 
Eqs. ()36p and ()37|h we can see the difference of the correlation properties of the two states' 
group delays. Considerable attention should be paid to that the correlation of the entangled 
state will increase, instead of decrease, when the propagation distance of L increases. Such 
effects should be able to be observed by means of the experimental methods shown in 




V. CONCLUSIONS 



We have demonstrated some properties of the "mode-entangled states" as the classical 
simulation of the quantum entangled states. These properties can be regarded as the proofs 
of the existence of the "mode-entangled states" . Then two experimental schemes to demon- 
strate these properties are suggested. One is based on the violation of the Bell inequality, the 
other on the correlation properties of the optical pulses' group delay in random waveguides. 
As far as we know, both of the two schemes can be carried out in current experimental 
conditions. We are looking forward to performing relevant experimental schemes. 
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Fig. 1: Electric field profiles for the optical transverse modes TE (the symmetric mode) 
and TEi (the antisymmetric mode). 

Fig. 2: The intensity variances of Y-splitter's two branches by changing the RI of the core 
layer (the length L = 1mm): (a) An = 0, (b) An = 0.0001, (c) An = 0.00021. 

Fig. 3: Experimental scheme. The input of control field is given at the mode superposi- 
tion ^ (|TE ) + |TEx)) and the input of target field is given at the mode |TE ) or 
|TEx). Then the output fields of the C-NOT gate are sent to spatially separated mode 
analyzers, each of which contains a Y-splitter and a variable phase control 6\ (9 2 ). 
The detected photocurrents are passively subtracted and monitored on a spectrum 
analyzer (SA) to check for correlations. 

Fig. 4: BPM simulation result for the scheme shown in Fig. 3. 
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Fig. 1: Electric field profiles for the optical transverse modes |TE ) (the 
symmetric mode) and iTEj (the antisymmetric mode). 



Fig. 2: The intensity variances of Y-splitter's two branches by changing the RI of 
the core layer (the length L-lmm ): (a) Arc = , (b) An = 0.0001 , (c) 
An = 0.00021. 




Fig. 3: Experimental scheme. The input of control field is given at the mode 
superposition -^-(|TE () ) + |TE 1 )) and the input of target field is given at the mode 

|TE ) or |TEj) . Then the output fields of the C-NOT gate are sent to two spatially 
separated mode analyzers, each of which contains a Y-splitter and a variable phase 
control l (0 2 ). The detected photocurrents are passively subtracted and monitored 
on a spectrum analyzer (SA) to check for correlations. 



Fig. 4: BPM simulation result for the experimental scheme shown in Fig. 3. 



